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Abstract 

The integrability properties of the field equation L^x = F{x)L'^ of a spherically 
symmetric shear-free fluid are investigated. A first integral, subject to an integra- 
bility condition on F{x), is found, giving a new class of solutions which contains the 
solutions of Stephani (1983) and Srivastava (1987) as special cases. The integrability 
condition on F{x) is reduced to a quadrature which is expressible in terms of elliptic 
integrals in general. There are three classes of solution and in general the solution of 
Lxx = F{x)L?' can only be written in parametric form. The case for which F = F{x) 
can be explicitly given corresponds to the solution of Stephani (1983). A Lie analysis 
of Lxx = F{x)L'^ is also performed. If a constant a vanishes, then the solutions of 
Kustaanheimo and Qvist (1948) and of this paper are regained. For a 7^ we reduce 
the problem to a simpler, autonomous equation. The applicability of the Painleve 
analysis is also briefly considered. 



1. Introduction 



Spherical symmetry is often assumed when seeking exact solutions to the Ein- 
stein field equations for a shear-free matter distribution. Spherically symmetric 
solutions may be used as isotropic cosmological models without the further restric- 
tion of homogeneity, or to model the interior of expanding or contracting spherical 
stars (Santos 1985, de Oliveira et al 1985). The shear-free condition substantially 
simplifies the field equations and not surprisingly most of the exact solutions known 
are not shearing (see Kramer et al 1980, Krasinski 1989). 

For a spherically symmetric, shear-free perfect fluid the field equations reduce 
essentially to two nonlinear partial differential equations. Various approaches have 
been followed in seeking solutions to these equations. McVittie (1933, 1967, 1984) 
assumed a functional form of the metric coefficients which enabled him to solve the 
field equations to give a wide class of solutions. The solutions found by Kustaanheimo 
and Qvist (1948), Chakravarty et al (1976), Wyman (1976) and Stephani (1983) 
depend on a suitable choice of a function of integration. A more recent approach 
by Hcrrera and Ponce dc Leon (1985), Dyer et al (1987), Sussman (1989), Maartens 
and Maharaj (1990) and Maharaj et al (1991) is to suppose that the spacetime is 
invariant under a conformal Killing vector, which leads to classes of solutions. 

One of the field equations for a spherically symmetric shear-free fiuid can be 
reduced to the partial differential equation L^x — F{x)L'^, where F{x) is an arbi- 
trary function resulting from an integration process. The purpose of this paper is 
to investigate integrability properties of this equation. In section 2 we give the field 
equations for a shear-free matter distribution. A condition on the function of inte- 
gration F{x) is found in section 3 that yields a first integral of L^x — F{x)L'^. We 
express this condition as a third order ordinary differential equation which in general 
has solutions in terms of elliptic integrals. The first integral found does not contain 
the solution of Kustaanheimo and Qvist (1948), but docs contain those of Stephani 
(1983) and Srivastava (1987) as special cases. This represents a new class (as far as 
we are aware) of solutions, generalising those of Stephani. In section 4 we perform a 
Lie analysis on this equation. If a constant a vanishes, then the solution is reduced 
to a quadrature. One class of the a = solutions is just the Kustaanheimo and Qvist 
solution. The other class is that found in section 3, thus giving a basis for the success 
of the ad hoc method used in section 3. For the case a 7^ we reduce the problem 
to a simpler, autonomous equation. In section 5 we consider the possibility of using 
the Painleve analysis and discuss results obtained for F{x) — x". One consequence 
of the results is to situate the known solutions and understand what is special about 
them. 



2. Field equations 



For a spherically symmetric, shear-free, perfect fluid we can introduce a comoving 
and isotropic coordinate system a;* = {t, r, 9, 0) such that the metric can be written 
as 
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Under these conditions the Einstein field equations reduce to (Kramer et al 1980) 

e'^ = A,e-^(*) (2) 
e^(A,, - - A,/r) = -F(r), (3) 

where f{t) and F{r) are arbitrary functions of integration. The energy density jj, 
and pressure p then assume the form 

/i = 3e^f -e-^^{2Xrr + \l + 4\/r) (4) 
p = \;^e-'^dt[e\Xl + 2Xr/r)-e'^+'^)]. (5) 

To find an exact solution of the field equations (2)-(3) the functions f(t) and F{r) 
need to be specified and equation (3) has to be solved for A. Then the dynamical 
quantities /x and p can be computed from (4) and (5). 
The transformation (Kustaanheimo and Qvist 1948) 

X = L{t, x) = e"^, F{x) = F/Ar^ 

reduces equation (3) to 

L,, = F{x)L\ (6) 

We will be concerned with this form of the field equation in subsequent sections. 
In spite of the tremendous interest generated by spherically symmetric shear-free 
perfect fluids the solution of the field equations in complete generality is known only 
for (Kustaanheimo and Qvist 1948) 

F{x)^{ax'' + bx + cy^^^ (7a) 

and (Stcphani 1983) 

F{x)^x-'^/' or (76) 

Incidentally subclasses of the McVittie (1984) class of metrics also correspond to 
these forms of F{x). We extend this to a new class of solutions in section 3. 

As pointed out in Kramer et al (1980), the invariance of (6) under coordinate 
transformations that preserve the form of (1), implies that F transforms as 

F(x) F(x) = x-^F{x-^) (8a) 

under 

X ^ X — x"^, L ^ L — x~^L (86) 
(Note the erroneous x^, instead of x~^, given on p 168 of Kramer et al.) This follows 



smce 



Thus a solution for a given F{x) is equivalent to a solution with F{x) replaced by 
x'^F{x~^). Note that (7a) is form invariant under (8), while the pair in (76) are 
equivalent. 

Note also that the solutions of (6) with F = {ax + 6)" are trivially related to the 
solutions with F = x". However, the two types of solution are not equivalent, since 
the metric is not invariant under x ^ ax + b (Stephani 1983, Srivastava 1987). 
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3. A first integral of L^x — F{x)L'^ 

Here we extend the technique of Srivastava (1987) to obtain a first integral of (6) 
without choosing an exphcit functional form for F{x). The first integral found is 
subject to an integral equation in F{x) which can also be expressed as a nonlinear 
third order ordinary differential equation. 

A number of choices of F{x) is known that reduces the order of the field equation 
(6) (Srivastava 1987, Stephani 1983). Rather than choose F{x) we seek a general 
condition on F{x) that reduces (6) to a first order differential equation. The integral 
of F{x)L'^ can be expressed as 

j F{x)L'^dx = FiL'^ -2 j FiLLJx, (9) 

where for convenience we have used the notation 

Fr = j F{x)dx. 

Integrating FiLL^ by parts and using (6) gives the result 

J FiLLxdx = FiiLLx - j FuLldx - j FFuL^dx 

On further evaluating the integrals on the right hand side of this equation and sub- 
stituting in (9) we obtain 

J F{x)LHx = FiL^ - 2FiiLLx + 2FiiiLl + 2{FFn)iL^ 

-2 \ j [2FFni + ?>{FFn)i\ L^L.d J (10) 

We observe that the integral on the right hand side of (10) can be evaluated if 
2FFiii + 3{FFii)i is a constant. Then equations (6) and (10) yield the result 

Mt) = -L. + FiL" - 2FnLL, + 2FniLl + 2 [(FFn)i - |i^o] (11) 

when we impose the condition 

2FFin + 3{FFn)i^Ko, (12) 

where "00 (^) is an arbitrary function of integration and Kq is a constant. Thus we 
have established that (11) is a first integral of the field equation (6) provided that 
condition (12) is satisfied. 

We transform the integral equation (12) into an ordinary differential equation 
which is easier to work with. Differentiating equation (12) gives 

2F,Fni + hFFn = Q (13) 

The form of (13) suggests the transformation 

:F = Fni. (14) 
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Substitution of (14) into (13) gives 



with solution 

T^xx = K,T-'''\ (15) 

where Ki is a constant of integration. Note that i^i = imphes F = by (15) and 
(14). We have estabhshed that the third order ordinary equation (15) with definition 
(14) is equivalent to the integrabihty condition (12). A solution of (15) gives T and 
then F{x) is found from (14). Repeated integration of (15) yields the result 

^K^ + Ks(^l + (/ J'-'^'dx^ ' - ^K, J^-'/^dx^\ (16) 

where K2,K3, are constants of integration. As this integration is nontrivial we 
provide details in the appendix. At this point it is convenient to let 



u 



^Jr-'/'dx 



{T-^f. (17) 



so that 

Ux 

Substitution of (16) into (17) gives the result 

x-xo^ [ — (18) 



{K^ + K^u + K2U^ - \Kiu 

where Xq is constant. Thus wc have reduced the differential equation (15) to the 
quadrature (18) which can be evaluated in terms of elliptic integrals in general (Grad- 
shteyn and Ryzhik 1994). In summary: the first integral (11), with F given by (18), 
represents a class of solutions of (6). 

It would be of interest to identify those cases where we can perform the inversion 
u = u{x) and explicitly write F = F{x): Once the integral in (18) is evaluated we 
invert to obtain u = u{x) and then (17) gives = u~'^l'^ . Finally the function F{x) 
is found from F — Txxx- This F then satisfies the integrabihty condition (12) for 
the existence of the first integral (11). 

Hence to find solutions F{x) in closed form, satisfying the integrabihty condition 
(12), we need to evaluate the integral (18). It is convenient to express our solution 
to (12) in the parametric form 

F{X^ J' j;xX 

Ux = J'-'/'^Wiu)]-' 

x-Xq = g{u), 
where for convenience we have set 

du 
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(u) = / , 

{Ki + K^u + K2u'^ - ^K^u^ f 
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Three cases arise in the above solution depending on the nature of the factors of the 
polynomial K4 + K^u + K2u'^ — \Kiu^: 

Case I: Three repeated factors 

If there are three repeated factors then we can write 

+ K^u + K2U^ - ^Kxu^ = (a + huf b ^ 0. 

In this case (18) can be integrated immediately in terms of elementary functions and 
we obtain 

»(«) = ~(a + bu)-y' 

Thus in this case it is possible to invert the integral (18) and to write u — u{x). 
After reparametrisation we can write 

F{x) = X-''/' 

which is the same as the case reported by Stephani (1983). We recall that this also 
leads to the a;"^'^/'' solution as noted in section 2. 

Case II: Two repeated factors 

With two repeated factors we can write 

+ K3U + K2u'^ - ^Kiu^ ^{a + bu){u + cf b 

The function g{u) in (18) may be easily integrated for two repeated factors so that 



1 56 156^ \ 1 

156^ r du 



^^""^ [ 2{a - bc)v? ^ 4(a - bcfu ^ 4(a - bcf ) y/a-bc + bu ^ 



[ 

J u\/a 



8{a - bc)^ J uVa -bc + bu' 

where the integral on the right side is expressible in terms of elementary functions. 

(The precise form depends upon the sign of a — be. See Gradshteyn and Ryzhik 
(1994) 2.246 and 2.248.9.) It is clear from the above form for g{u) that it will not 
be possible to perform the inversion to obtain u = u{x). (This is possible only when 
b — which is not permitted.) Consequently the solution can only be given para- 
metrically in this case. 

Case III: No repeated factors 

If there are no repeated factors, we write 

+ K3U + K2u'^ - ^Kiu^ = d{a -u){b- u){c - u) d^O. 
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In this case the solution of (18) is given in terms of eUiptic integrals. We obtain 

, , , o/n 2\c(a — c) + b(a — b) — u(2a — c — 6)1 
g{u) = d'"^'^ ^ / ^ ^' + 



a — b)(a — c)(b — c)^^ (a — u){b — u){c — u) 

^-3/2 1 ^_ X 



(a-6)2(6-c)V(«-c 



^3 



X 



(6 - c)(a + 6 - 2c)F(q;,p) - 2{& + + b-ab-ac- bc)E{a,p) 



where we have set 



la — b 



a — arcsmW p — 

\ a — u 



In the above F{a,p) is the elliptic integral of the first kind and E{a,p) is the elhptic 
integral of the second kind. By inspection of g{u) we observe that it is not possible to 
invert and obtain u = u{x). Therefore, as in Case II, the solution may be expressed 
only in parametric form when there are no repeated factors. 

Wc have established that for our class of first integrals the form of F may be given 
in parametric form in general. It is only in Case I that inversion is possible and we 
can write F — F{x) exphcitly. This case corresponds to F = x~^^^'^ which was first 
identified by Stephani (1983). 

It is interesting to observe that the solution of Kustaanheimo and Qvist (1948) 
is not contained in the first integral (11). To prove this wc simply substitute (7a) 
into (12) and see that (12) is not satisfied. Thus the first integral (11) cannot admit 
the Kustaanheimo and Qvist (1948) solution with (7a). The first integral (11), with 
F given by the quadrature (18), represents a new class of spherically symmetric 
solutions. Most papers have dealt with special cases of (7a) (Kramer et al 1980). 

Note that the first integral (11) was obtained without specifying the function 
F{x). However, in our solution F{x) is constrained by the integral equation (12). 
Thus any F{x) satisfying (12) will yield a first integral of the form (11). The solutions 
of Srivastava (1987) and Stephani (1983) can be regained as special cases of (11). 
With the choice F{x) = {ax + 6)" where n — —15/7 we find that (12) is satisfied 
with Kq = 0. Then (11) becomes 

Mt) = imMt) 

= -6L^ - (21/ Aa)(ax + b)~^'^ - (3/2)(7/a)^(ax + b)'^'^ LL^ 

+{l/A){l/af{ax + bf'Ll - (l/6)(7/a)'(ax + b)-^" L\ (19) 

The first integral (19) was reported by Srivastava (1987). Now with a = 1 and 6 = 
equation (19) reduces to 

99o(t) = (l/6)Vo(t) 

= - (21/4)x-«/^L2 - (3/2)72a;-i/^LL, + (1/4)7V/^L2 

-{l/Q)tx-^/'L\ (20) 

The first integral (20), where F{x) is given by (8), was found by Stephani (1983). 
Note that in (19) and (20) the functions 0o(^) and ^poit) have been introduced to 
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ease comparison with the corresponding expressions given by Srivastava (1987) and 
Stephani (1983) respectively. Hence the first integral (11) represents a class of solu- 
tions which may be interpreted as a generalisation of the Stephani/Srivastava solu- 
tion. In fact, the F of the form [ax + 6)" only satisfies (12) if n = —15/7 (which can 
be mapped to the case n — —20/7 by (8)). This is a special solution of (12). The 
general solution, containing the generalisation of Stephani and Srivastava, is given 
by (18). The case n = —15/7 is characterised by the fact that it is the only case 
for which we can write F = F[x) explicitly; all other solutions may only be given 
parametrically. 

4. Lie Analysis 

In section 3 we obtained the first integral (11) of (6) in an ad hoc fashion. Here 
we consider a widely applied technique to reduce (6) to a first order differential 
equation. A systematic method to determine whether a second order ordinary dif- 
ferential equation can be solved by quadratures is that of Lie (1912). On performing 
the Lie analysis (Leach 1981, Leach et al 1992) we find that (6) is invariant under 
the infinitesimal transformation 



provided that 



F[x) = M^-^/^exp(|/^) (22) 
A_ = 2MC^-^/^exp(^||^j (23) 



B{x) = {A,-a)/2 
C{x) = Co + Cix, 

where a, M, Co and Ci are constants. A solution of (23) gives F{x) by (22) which 
reduces (6) to a first order differential equation. Note that the existence of the 
symmetry (21) permits a reduction to first order. If two symmetries G are known, 
the second order equation (6) may be reduced to quadratures. 

For convenience we divide the integration of (23) into a number of cases: 

Case r.a^O C = 

This is the simplest case. Then (22) yields 

F{x) = M{ax^ + hx + c)'^''^ 
which is equivalent to the Kustaanheimo and Qvist (1948) class of solutions. 

Case II: a = Co ^ = Ci 
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In this case (23) reduces to the differential equation 

where Ki = 2MCq. This is the same as equation (15) considered in section 3 and 
the results of that section become applicable. Thus we have estabhshed that the first 
integral (11) is a special class of solutions admitted by the general Lie method. 

Case III: a = Ci 7^ 

We now consider the case Ci 7^ 0. We begin by observing that (23) is simplified 
by a change of the independent variable x which transforms away the linear term 
C{x) — Cq-\- Cix. If we introduce the canonical variables (Ci 7^ 0) 

Y = A/C\ X = -l/C, 

(23) is transformed to 

Y"' = MiY-^'"" (24) 

where Mi = 2MC]"^, and we have used the notation ' := d/dX. We provide an 
alternate solution to (24) than that presented in section 3. Equation (24) can be 
reduced to an equation involving only Y' if we introduce the variable 

s := J Y-^'^dX (25) 

Then (24) and (25) imply 

YY's' = -M3 - 2M2S + |MlS^ 

where M2 and M3 are constants. On multiplying this equation by s' and integrating 
we obtain 

y^s'^ = M4 + M3S + M2S^ - |MlS^ (26) 
where M4 is constant. Equation (26) is reduced to the quadrature 

x-xi^ [ — (27) 

(M4 + M3S + M2S2 - |Mis3)'^/2 

where Xi is a constant of integration. As mentioned in section 3 integrals of the 
form (27) can be evaluated in terms of elliptic integrals in general in which case the 
solution will be given parametrically. However, in some cases after having solved 
(27) we then invert to obtain s{X) and finally we have A — C^Y — C^{s')^'^. 

Case IV:ay^O C = 

It remains to consider the highly comphcated case of a 7^ 0. Apart from the 
subcase C = 0, we have not been able to integrate (23) which is hardly surprising 
considering its highly nonlinear nature. The special case is a ^ = C. This appears 
to give a generalisation of Kustaanheimo and Qvist since (22) and (23) imply 

A — ax^ + bx + c 



F = MA-^/'' exp (^J ^dxy 
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Now let 



Then (6) is autonomised: 



y" - ay' +(^ + ac- ^bA y + uy^ = o. (28) 



Now, if A = 6 — Aac > and a = ±5v A, then (28) can be transformed to 
where 

r] — e' 

Then (29) is easily integrated to yield 

drj 



2aX/5y 



[C- (2M/3A)r/3]V2- 
This seems to be a new solution, but 

A > ^ ^4 = ax"^ + bx + c — a{x — fj,){x — v) 
with jji — y = y/A/a (t^ 0). Then 

dx 1 , /x — u 
m 



A x/A \x - u 



implies 

[ {Ka-'/^){x - z/)-5 a = -5y/A 

which remarkably is a special case of Kustaanheimo and Qvist! Thus even with 
Q! 7^ = C we do not obtain any new solutions. 

We note that if instead of assuming that a has the value that enables us to 
transform (28) to (29) other functions F can be obtained. For example if we put 
a = 6 = in ^4, 

F{x) = Mc-^/2e"^/(2^) 
which is equivalent to the Stephani (1983) form 

F{x) = 
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after reparametrisation. When A — a{x — P){x — Q), 

F{x) = Ma^^'"^ {x ~ P)2^^W=Q)~^ (x - Qy^^(^-^ 
which reduces to to the Stephani (1983) form 

F{x) = (x + q;)"(x + /3)-"-^ 
on the substitution 2a{p-Q) ^^+2- If^ = a(3^ — -P)^ then 

F{x) = {x - Py^ e-^^ 

which via (8a) is related to F{x) — e^. 

Case V: a^O Ci ^ 

This is the most difficult case to analyse. When Ci ^ 0, we are in a position to ex- 
press (23) as a simpler autonomous equation after a change of variable. Introduction 
of the variable 

a f dx^ 



^ = exp-y- (30) 



transforms (23) to the differential equation 



- 



2{5/af/^MCA-^/^. (31) 



A J 

For Ci 7^ without any loss of generality we can write (31) as 

X) = ^^'^'''^ (32) 

•^X' XXX 

where we have set Co — and 2{h/af^'^Ci — 1. Essentially we have rescaled the x 
variable by x — > Px + 7. Expansion of the left hand side of (32) results in 

XXX ^AAxAxxAxxx ^^x'^xx ^•^•^xx •^■^x ^ — 0. (33) 

The fourth order equation (33) is reduced to autonomous form if we introduce the 
new independent variable 

1 

X 

With this transformation (33) assumes the autonomous form 

The order of (34) is reduced by the substitution 

A = e\ Ar = AZ^e-^' 

to 

ZZ,,, - Z,Z,, - ^ZZ,, + bZl + = 0. (35) 

Thus we have expressed equation (23) with a 7^ in the simpler form (35). A 
solution of (35) gives A and A is then found from (30). 
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5. Discussion 



We have considered the integrabihty of the equation L^x = F{x)L'^. A first 
integral was found in section 3 that contains the special cases of Stephani (1983) 
and Srivastava (1987) studied previously. The solution is given parametrically and 
there are three classes of solution. In one case wc can perform the inversion and 
write F = F{x) explicitly. This corresponds to F = x^^^^^ obtained by Stephani. In 
section 4 the Lie analysis of this equation resulted in the equation (23). For a = 
equation (23) has solutions in terms of elliptic integrals. In the more complicated 
case of q; 7^ we reduced (23) to the simpler autonomous equation (35). 

Considering the significance of spherically symmetric models in general relativity 
the equation L^x = F{x)LP' deserves further attention. Another approach in the 
investigation of the integrabihty of this equation is the Painleve analysis. In the 
Painleve analysis one performs a singularity analysis of a differential equation. If 
the only movable singularities are poles, the equation is said to possess the Painleve 
property and it is conjectured that such equations can be solved analytically (Ramani 
et al 1989). 

Omitting the details we illustrate this approach with an example . With F[x) = 
x"' Stephani (1983) obtained an ansatz that reduced 

Lxx = x'^L^ (36) 

to a first order equation. We performed the Painleve analysis on equation (36) 
after transforming it to autonomous form by the transformation which converts the 
appropriate form of G in (21) to d/dX. The autonomous form of (36) is 

Y" - (2n + 5)Y'+{n + 2) (n + 3)Y -Y^ = (37) 

and satisfies the Painleve property for the cases n = —5/2, 0, —5, —15/7 and —20/7. 
These are precisely the cases of known solutions, i.e. Kustannheimo and Qvist (1948) 
(n = -5/2, 0, -5) and Stephani (n = -15/7, -20/7). 
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Appendix 

In this appendix we indicate how the third order equation (15) may be integrated 
to yield the solution (16). The nonlinear equation 

■r ly' -77—5/2 

XXX -IV 1^ 

may be written as 

{^^.^)x - \ (•^.'), = ^1-^"'/' 

This gives 

where —IK^ is a constant. This second order equation may be written in the form 
On integration we obtain 



2 



where —Kj, is a constant. It is convenient to rewrite this first order equation in the 
form 



A final integration generates the solution 



2 

X 



where is a constant. 
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